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Abstract. The similarity between the mass and spatial distributions of pre- 
stellar gas cores in star-forming clouds and young stars in clusters provides 
strong circumstantial evidence that these gas cores are the direct progenitors of 
individual stars. Here I describe a physical model for the evolution of massive 
cores into stars, starting with the intial phases of collapse and fragmentation, 
through disk formation and fragmentation, the later phases of stellar feedback, 
and finally interaction of the newly formed stars with their environments. This 
model shows that a direct mapping from cores to stars is the natural physical 
outcome of massive core evolution, and thereby allows us to explain many of the 
properties of young star clusters as direct imprints of their gas-phase progenitors. 



1. Introduction 



Massive stars form in regions of extremely high column density, hidden behind 
hundreds of magnitudes of visual extinction. As radio and submillimeter inter- 
ferometers have matured over the past decade, they have revealed the properties 
of these regions with ever higher detail, to the point where today it is not entirely 
unreasonable to speak of observationally determined "initial conditions" for the 
problem of star formation. One of the most striking results of this exploration 
has been the extent to which the properties of young star clusters are directly 
mirrored in the conditions found in pre-stellar molecular gas. 

As the most basic level, young clusters and the molecular clumps from 
which they form have similar bulk properties, such as column density 1 g 
cm~^ ), size (~ 1 pc), and velocity dispersion (a few km s~^) (jMcKee fc Tan 
20031 ) . More interestingly, the dense cores within these clumps also mirror the 
properties of stars. Cores are bound, centrally-condensed objects with charac- 
teristi c sizes ~ 0.1 pc or sinaller and masses co mparable to those of individual 
stars ( Sridharan et al.ll2005l : iBeuther et al.ll20()5l ). Observations in many regions 
with a variety of techniques find that the core and star mass functions have 
the same shape, differing only in that cores are a factor of 2 — 4 more massive 
(e^ lMotte et al.lll998l : iJohnstone eralll200ll : iReid fc WilsOT]l2006l : lAlves et al 



20071 . and J. Alves, this volume). Moreover, cores appear to be mass-segregated: 
cores with masses greater than a few Mq are found only in the centers of their 
parent clumps, but the core m a ss function is other wise independent of position 
( Elmegreen &: Krakowski 2001 : Stanke et al. 20061 ). This is remarkably similar 
to the pattern in young clusters, where there is no segregation for stars smaller 
than a few Mfn , but more massive stars are almost ex clusively in cluster centers 
(iHillenbrand fc Hartmaiinlll998l : iHuff Stahlei1l2006l ^. 
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The coincidence in both the mass and spatial distributions of cores and 
stars makes a strong circumstantial case that young stars' properties might be 
direct imprints of core properties. (The somewhat higher core masses are to be 
expecte d, since outflows should p revent ~ 50% of a core's mass from reaching 
a star - iMatzner McKed 12000'.) However, we cannot make such an inference 
without a theoretical understanding of how cores might evolve into stars. The 
goal of this paper is to summarize work in the last few years that has focused 
on the most problematic part of this process, understanding the evolution of 
massive cores. In the following sections I sketch a model for the evolution of 
these objects, beginning with observed core properties and using numeric and 
analytic arguments to understand how they collapse into stars. 



2. Initial Collapse and Framgentation 

The first phase of evolution for a massive core begins when it starts to collapse 
but has not yet formed any stars. One might initally suspect that there is no 
plausible way for a massive core to collapse to a single star or a small-multiple 
system, since the Jeans mass in cold molecular clumps is only ~ Mq. If massive 
cores do truly fragment into Jeans mass-sized objects once their collapse begins, 
then there can be no direct mapping from cores to stars. Such behavior is exactly 
what purely hydrodynamic numerical simulations find: as objects collapse and 
their density rises, the Jeans mass falls, and the objects break into smaller 
and smaller pieces that always have masses comparable to the Jeans mass at 
their current density. The fragmentation process ceases only when the assumed 
equation of state stiffens, so massive co res generate numerous small stars (e.g. 
Bate fc Bon^l2005l : iDobbs et al.ll2005l ). 



However, while a purely hydrodynamic approach to fragmentation is analyt- 
ically very simple and numerically very cheap, it neglects th e importaii t effect 
of radiation feedback from embedded, forming protostars. KrumhoI3 ( 20061 ) 



points out that, even before embedded stars begin nuclear burning, just the 
gravitational energy released as gas accretes onto them can significantly heat 
the surrounding gas, r aising the Jeans mass and suppressing fragmentation. 
Krumholz et al. ( 2007d ) follow up this point by simulating the collapse of tur- 



bulent, massive cores whose initial masses, sizes, central concentrations, and 
lev els of turbulence are c hose n to match those of observed cores (as seen e.g. 
by ISridharan et al. 2005 and Beuther et al. 20051 ). The simulations combine 



a protostellar evolution model with a new adaptive mesh refine ment gravity- 
radiation- hydrodynamics algorithm ( Krumholz et al.]|2004 . 2007al ) to model ac- 



curately the effects of radiative heating. The simulations show that radiative 
heating strongly suppresses fragmentation of massive cores, allowing the great 
majority of the mass in a core to collapse into one or two stars. In contrast, a 
control run omitting radiative heating qualitatively reproduces the earlier hy- 
drodynamic result that massive cores collapse into dozens of small fragments. 
Figure [T] illustrates the difference made by the inclusion of radiative heating. 

The most important point to take from this is that, when more detailed 
physics than simple hydodynamics is included, simulations and analytic argu- 
ments both show that the observed massive cores are unlikely to fragment into 
many pieces. Thus, a model in which there is a direct mapping from the masses 
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Figure 1. Column density in two simulations of the collapse and frag- 
mentation of a 100 Mm core. The left and right panels show two runs from 
iKrumholz erall (|2007cL runs lOOA and lOOISO) at 20 kyr of evolution; one 
uses radiative transfer and one does not, but otherwise the simulations have 
identical initial conditions and resolution. The plus signs indicate the posi- 
tions of stars. Note the significantly greater number of stars, with several 
more condensations on the verge of collapse to stars, in the non-radiative run. 



and positions of massive cores to those of massive stars passes its first test: the 
cores will collapse largely monolithically, rather than fragmenting to small stars. 
A secondary point is that we must be careful about drawing conclusions based 
on simulations with very simple physics. Models that do not include radiative 
transfer produce qualitatively different results from those that do. 



3. Disk and Binary Formation 



Since collapsing turbulent cores have non-zero angular momentum, they nat- 
urally form protostellar disks. Observationally, these are a potential signpost 
of the star fo rmation process. Both th e simulations described above and ana- 
lytic models ( Kratter Sz Matzneil l2006l ) find that the disks formed by massive 
cores are likely to be strongly gravitationally unstable. This instability causes 
the disks to develop large-amplitude m = 1 spiral modes, and potentially even 
fragment to form companions to the primary star (al though a majority of the 
mass still goes into the primary, not the fragments). IKrumholz et al.| ( 2007bl ) 
shows that the strong m = 1 spiral structure present in such unstable disks 
should be observable with next-generation telescopes such as ALMA and the 
EVLA, and that a systematic offset between the disk's "zero" velocity and the 
central star's velocity produced by the instability might also be detectable. The 
detection of a disk with these signatures inside a massive core would be strong 
evidence in favor of the model that massive cores collapse to form massive stars. 
Figure [2] shows a simulated ALMA obse rvation of such a disk, computed using 
the technique of Krumholz et al. ( 2007bl ). 
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Figure 2. The velocity-integrated brightness temperature (left panel) and 
centroid velocity {right panel) in a simulated 0." 1 resolution ALMA observa- 
tion of a massive protostellar disk 500 pc away in the CH3CN 220.7472 GHz 
line. The disk is from the simulation illustrated in the left panel of Figure [T] 
at an evolution time of 27 kyr, when the central star is 8.3 Mq in mass. The 
white contour in the right panel corresponds to —5 km s~^, and the black 
contour shows km s~^, measured relative to the velocity of the central star. 
The disk is systematically offset to negative velocities relative to the star, so 
while there is a large region with velocity < —5 km s~^, there are no pixels 
with velocities > 5 km s^^. Black pixels correspond to locations where ALMA 
would not detect the line at > 3cr confidence. 



Disks around massive stars are also important for their role in fo rming 
comp anions to massive stars. Most massive stars have close companions (|Lada 



20061 ) ■ and gravitational instability is massive protostellar disks provides a nat- 



ural explanation for this because even radiatively heated massive disks suffer 
some fragmentation (though vastly less than if radiation is omitted). These 
fragm ents initially form with m asses ~ M p) at distances > 100 AU from the pri- 
mary ( Kratter &; Matzneil 12006, : .Krumholz et al.ll2007d ). but they subsequently 
migrate inward to separations < 10 AU as the disk accretes. The final separa- 
tions of these disk-formed companions from the primary, and whether some of 
them merge with it, has not yet been determined. 

An interesting fate awaits migrating stars that get close to the primary 
but do not merge with it. Massive protostars go through a phase of deu- 
terium shell burning, during w hich their radii swell to tenths of an AU in size. 
Krumholz Sz Thompson ( 20071 ) point out that this may lead a primary to over- 



flow its Roche lobe and transfer mass onto close companions. Such mass transfer 
is almost always unstable, terminating only once the mass ratio of the system 
reaches unity. This provides a natural mechanism for the origin of the hereto- 
fore unexplained massive "twins" , binarie s consisting of two massive st ars with 
a mass ratio of almost exactly unity (e.g. Pinsonneault Sz Stanek 20061 ). 
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4. Radiation Pressure Feedback 



A direct core to star mapping is possible only if most of a core's mass be able 
to accrete onto the massive star forming within it. However, spherically sym- 
metric calculations indicate that the huge radiation output output of massive 
stars should exert a force stronger than gravity on dust grains suspended in 
the gas around them for stars of ~ 20 MfD or lar ger (jLarson Sz Starrfieldlll971 



lYorke Kru"e rei"l977l: IWolfire k Cassinellilll987l ). How massive stars can form 
despite this radiation barrier is a classic problem in astrophysics. Fortunately, 
rotation and the formation of a disk can significantly mitigate this effect, because 
gas in a disk self-shields agjainst the radiation du e to high optical depths ( Nakanol 
19891 : iNakano et al.|[T995l : Ijiiina k AdamslfTOQ^ ). while at the same time a disk 



collimates the radiation field and beams it away preferentially in the polar di- 
rection, thereby reducing t he radation force felt by gas in the equatorial plane 
(lYorke k Sonnhalted [200I 1 . However, even with these effects it is not entirely 
clear that stars can grow to arbitrary masses by accretion. 

Two additional effects may help. First, radiation hydrodynamic simula- 
tions show that the first effect of radiation pressure is that massive stars blow 
radiation bubbles abov e and below accretion disks ( Yorke k Sonnhalteil |200 
Krumholz et al.ll2005al ). However, in three dimensions these bubbles do not halt 



accretion, because gas that reaches the bubble wall flows along the wall onto the 
accretion disk. Bubbles may also collapse due to Rayleigh- Taylor instability, 
allowing accretion to continue through optically thi ck channels while radia, tion 
escapes through optically thin regions around them ( Krumholz et al.]|2005al ). If 
magnetic fields are present and sufficiently str ong, this effect will be enhanced by 
photon bubble instability ( Turner et al. 20071 ). which arranges the gas into dense 
lumps separated by low-density gaps through which radiation leaks, effectively 
reducing the radiation pressure force experienced by the bulk of the gas. 

Second, protostellar outflows provide a third escape valve for radiation. 
Massive protostars appear to generate hydromagnetic outflows just like low mass 
stars, with the difference that for massive stars the outflow cavities are largely 
dust-free because the base of the outflow is close enough to the star for the 
dust within it to have been destroyed by sublimation. Such outflow cavities 
therefore present optically thin channels through which radiation can leak out 
of the optically thick cores. Radiative transfer calculations show that this can 
lead to order-of-magnitude reductions in the radiation force on accreting gas 
near the equatorial plane, a gain allowing accretion to continue where it might 



otherwise have been halted ( Krumholz et al.ll2005bl ). 



While a definitive numerical simulation including the effects of radiation 
forces, magnetic fields, and protostellar outflows in three dimensions has not 
yet been done, and is probably at best barely within the capabilities of present- 
day supercomputers, it seems clear that each of these effects will help massive 
stars form by accretion. Thus, we can tentatively say that there is no barrier to 
most of the gas in a protostellar core accreting onto a massive star. Conversely, 
however, simulations of massive star formation that do not include radiation 
force effects are simply ignoring this problem entirely. Preliminary simulations 
indicate that models of star formation that depend on Bondi-Hoyle accre tion 
are likely to fail once radiation pressure is included ( Edgar k Clarkd |2004| ) . so 
models of this sort must be viewed with caution. 
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5. Competitive Accretion 

Thus far we have seen that massive protostellar cores will not fragment strongly, 
what fragmentation they do show is consistent with the observed multiplicity 
properties of massive stars, and that radiation pressure will not prevent most 
of the mass in a core from accreting onto a star. There remains, however, one 
more way in which a direct mapping from cores to stars could fail: if, once 
stars accrete their parent cores, they were subsequently to accrete a great deal 
more mass, then there would be no direct relationship between core and stellar 
masses. This process of accretion onto stars from gas that was not origina lly 
part of a bound core is known as competitive accretion ( Bonnell et al. 2001al lbl). 



Purely hydrodynamic simulations of star cluster formation show that this 
process is the dominant mechanism by which stars gain mass. In effect, all stars 
are born at masses of order the Jeans mass, but some of them fall to the center 
of the collapsing gas cloud, and the deep potential well then channels gas to 
them. They subsequen tly accrete this gas and grow in ma.ss, pr oducing a full 
range of initial masses (|Bonnell et al.ll2004l : [Bonnell Batd l2006l l . 



However, direct observational estimates of the rate of competitive accretion 
generally find that it is too sm al l to make a sign i ficant contribution to final 
stellar masses (|Andre et al.ll2007[ ). iKrumholz etal] (|2005d . l2006l ^ point out that 



competitive accretion is possibly in simulations only because the simulated gas 
clumps are in the process of global collapse, which creates deep potential wells 
within which the gas is dense and non-turbulent and Bondi-Hoyle accretion is 
rapid. These deep, dense, quiescent gas wells have not been observed, however. 
This is probably because clumps are not in a state of global collapse. Such a 
collapse necessarily converts order unity of the m ass in a gas clump in to stars and 
ends star formation in 1 — 2 free-fall times (e.g. i Bonnell et al.|[2004 )- However, 
there is strong evidence that the star formation process cannot be anywhere 
near that fast . The gas clump from which the ONC formed likely had a density 
10^ cm~^ ( Elmegreen I2OOOI '). implying a free-fah time of 0.1 - 0.2 Myr, but 



the estimated ages of the stars in Orion point to a formation process lasting 
1 — 3 Myr, implying a minimum formation time scale of 5 free-fall times even 
if one assumes the fastest pla usible formation time, with ~ 15 free-fall times 



being more likely (|Tan et al.l i2006). Furthermore, the total galactic mass of 
infrared dark clouds, which have densities ~ 10'^ cm~^, is such that, if they were 
collapsing to form stars on a free-fall time scale, the total galactic star formation 
rate would have to be ~ 100 times higher than its observed value. The same 
conclusion holds for dense gas clumps traced by HCN (1 — > 0) emission, which 



have densities ~ 10^ — 10^ cm ^ ( Krumholz Sz Tajil l2007'). This implies that, even 



in regions with mean densities of 10^ cm~^ star formation cannot be anywhere 
near as rapid as would be required for competitive accretion to occur. 

The final piece of evidence against c ompetitive accretion co mes from simula- 
tions t hat include more detailed physics. iLi &: Nakamura ( 20061 ) and lNakamura &: Li 



( 2OO7I ) simulate the formation of a star cluster using magnetohydrodynamics 



rather than simple hydrodynamics, and including the effects of outflows driven 
by the protostars forming within the cluster. They find that the protostellar 
outflows drive turbulent motions, preventing global collapse and ensuring that 
conditions are too turbulent for competitive accretion processes to alter the 
masses of protostars significantly after they have consumed their parent cores. 



Massive Cores to Massive Stars 



7 



Moreover, these simulations, unlike ones where competitive accretion occurs, 
produce star formation rates in good agreement with observations. We can 
therefore tentatively conclude that the core accretion hypothesis passes this fi- 
nal test: once stars have accreted their parent cores, they will not gain much 
additional mass from the gas to which they were not bound at birth. A direct 
core to star mapping will survive. 

6. Summary and Conclusion 

For the past decade observations have increasingly pointed to an intimate link 
between young stars and the dense stellar-mass gas clouds known as pre-stellar 
cores. Cores and stars have very similar mass and spatial distributions, so it is 
tempting to explain the properties of young star clusters as imprinted at birth. 
However, this hypothesis requires that cores map directly onto stars. Here I 
provide a physical model for such a mapping. When massive cores first collapse, 
they do not fragment strongly because the first stars that form within them heat 
the gas and suppress fragmentation. As a result, massive cores collapse to only 
a few stars. Most of the fragmentation that does happen occurs in unstable 
self-gravitating disks, which should be directly observable. Disk fragmentation 
ensures that massive stars will essentially always have companions, and some 
of these companions will turn into twins of the primary star. As the massive 
star grows, it will begin to generate a huge radiation force opposing accretion. 
However, a combination of instabilities and the leakage of radiation through 
protostellar outflow cavities allows the radiation to escape and accretion to con- 
tinue unimpeded to high masses. Finally, once stars have accreted their parent 
cores, they will be unable to gain additional mass that was not originally part 
of the core. Thus, the hypothesis that stars come directly from cores, with a 
one-to-one mapping of core mass to star mass, is in agreement with a physical 
model for how massive cores evolve. The best explanation for many properties 
of young star clusters appears to be that they are set at birth, when the cluster 
is still a dark cloud. 
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